Longitudinal boost-invariance of charge balance function 
in hadron-hadron and nucleus-nucleus collisions 
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Using Monte Carlo generators of the PYTHIA model for hadron-hadron collisions and a multi- 
phase transport (AMPT) model for nucleus-nucleus collisions, the longitudinal boost-invariance of 
charge balance function and its transverse momentum dependence are carefully studied. It shows 
that the charge balance function is boost-invariant in both p+p and Au+Au collisions in these two 
models, consistent with experimental data. The balance function properly scaled by the width of the 
pseudorapidity window is independent of the position or the size of the window and is corresponding 
to the balance function of the whole pseudorapidity range. This longitudinal property of balance 
function also holds for particles in small transverse momentum ranges in the PYTHIA and the 
AMPT default models, but is violated in the AMPT with string melting. The physical origin of the 
results are discussed. 

PACS numbers: 25.75.Gz,25.75.Ld 



I. INTRODUCTION 

Charge balance function (BF) has been widely used 
as an effective exploring for the hadronization scheme 
in hadron-hadron collisions at the ISR energies [l| and 
e"*" -|- e~ annihilations at PETRA energies 01- Re- 
cently, the charge BF gains special attentions in clocking 
hadronization at relativistic heavy- ion collisions. A nar- 
rowing of the BF is suggested as a signature for delayed 
hadronization [1, 

The dependence of the BF on centrality and sys- 
tem size has been reported by several relativistic heavy- 
ion experiments [1, 0|. However, most of the current 
heavy-ion experiments are limited by the pseudorapidity 
range 0, S 0] i it is impossible to quantitatively compare 
the results from the experiments with the coverage at 
different pseudorapidity ranges. The dependence of the 
BF on the pseudorapidity window is essential for under- 
standing the physics of the BF d, H, 1, [1] , and has been 
carefully studied by NA22 and STAR experiments for 
hadron-hadron and relativistic heavy ion collisions, re- 
spectively. 

The NA22 experiment has full Att acceptance and ex- 
cellent momentum resolution Q. It is found in the ex- 
periment that the BF in 7r+p and K+p Collisions at 22 
GeV is invariant under longitudinal boost over the whole 
rapidity range of produced particles, in spite of the non- 
boost-invariance of the single-particle density. Moreover, 
the BF of the whole rapidity range can be deduced from 
the BF properly scaled by the width of rapidity win- 
dows 0. 

The STAR experiment covers a finite but relative wide 
pseudorapidity range. The scaling property of the BF 
in Au -f Au collisions at 200 GeV is further observed in 
the experiment [l^l ■ This scaling property of the balance 
function is also found in different px ranges of final state 
particles. 

These results from both hadron-hadron and nuclear 
collisions indicate that charge balance of produced par- 



ticles in strong interactions is boost-invariance in longi- 
tudinal phase-space, in contrary with the single particle 
density. 

Therefore, it is interesting to see if those properties are 
taken into account in the models which are successfully 
described hadron-hadron and nuclear collisions, and how 
they associate with the mechanisms of particle produc- 
tion in the models. 



II. CHARGE BALANCE FUNCTION AND 
IMPLEMENT MODELS 

Charge balance function measures how the conserved 
electric charge compensate in the phase space, i.e., how 
the surrounding net charge are rearranged if the charge of 
a selected point changes In high energy collisions, the 
production of charged particles are constrained by charge 
balance in the phase space. The BF therefore provides a 
direct access to collision dynamics. 

The BF has been originally defined in terms of a com- 
bination of four kinds of charge-related conditional den- 
sities in pseudorapidity [l| 
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where UabiSr]) is the total number of pairs of opposite 
charged particles with pseudorapidity difference Srj in the 
pseudorapidity window ?7w and n_ are the number of 
positively and negatively charged particles in the window 
rjvt, respectively. (• • •) is the average over the whole event 
sample. 

From the findings of the BF at NA22 and STAR 
experiments [lo| . the BF is boost-invariant in the whole 
rapidity range in hadron-hadron collisions and may be in 
nuclear collisions as well. In the case, the properly scaled 
BF is corresponding to the BF of the whole pseudorapid- 
ity range and is deduced by 
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where |77w| is the width of pseudorapidity window. 

The PYTHIA 5.720 ^ is weh set up for p+p colli- 
sions. It is a standard Monte Carlo generator with string 
fragmentation as hadronization scheme. Two versions 



of a multi-phase transport (AMPT) model [IJI are used 
to study Au+Au collisions. One is the AMPT default 
and the other one is the AMPT with string melting. In 
both versions, the initial conditions are obtained from 
the HIJING model, and then the scattering among par- 
tons are given by ZPC. In the AMPT default model, the 
partons recombined with their parent strings when they 
stop interacting, and the resulting strings are converted 
to hadrons using the Lund string fragmentation model, 
while in the AMPT model with string melting, quark coa- 
lescence is used in combining partons into hadrons. The 
dynamics of the hadronic matter is described by ART 
model. 

It is commonly believed that in relativistic heavy ion 
collisions, the charge ordering during the string fragmen- 
tation in elementary collisions is no longer valid, and it 
should be replace d by the quark-coalescence mechanism 
in hadronization [l3| . So it is interesting to see whether 
the boost-invariance of the BF is sensitive to the mecha- 
nisms of hadronization. 

In this paper, we firstly study the boost-invariance of 
the BF for p+p collisions at = 22 GeV and ^/s = 200 
GeV using the PYTHIA, and for Au-|-Au collisions at 
^/s = 200 GeV using two versions of the AMPT. The 
transverse momentum dependence of longitudinal scal- 
ing property of the BF is then examined in the models. 
The obtained results are compared with corresponding 
experimental data and discussed. 



III. BOOST-INVARIANCE AND 
LONGITUDINAL SCALING OF THE BF 

In order to demonstrate directly whether the BF is in- 
variant under a longitudinal Lorentz transformation over 
the whole rapidity in hadron-hadron collisions, we choose 
four equal size (|/7w| = 3) pseudorapidity windows lo- 
cating at different positions (—3,0), (—2,1), (—1,2) and 
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FIG. 1: Upper panel: the B{5rj\ri„) in four pseudorapidity 
windows with equal size |77w| = 3 at the different positions 
for p+p collisions at (a) ^/s = 22 GeV and (b) -^s = 200 
GeV by PYTHIA model. Lower panel: the scaled balance 
function, Ba{Srj), deduced from the directly measured BF at 
six different sizes and positions of pseudorapidity windows for 
p+p collisions at (c) ^ = 22 GeV and (d) = 200 GeV 
by PYTHIA model. The solid down triangle is the BF of the 
whole rj range. 



(0, 3). The results for p+p collisions at -y/s = 22 GeV and 
-y/s ~ 200 GeV are shown in Fig. 1(a) and (b) respec- 
tively. The statistic errors arc smaller than the mark- 
ers. It is clear that the BF measured in four windows 
are approximately identical to each other at two inci- 
dent energies. This indicates that the charge compensa- 
tion is essentially the same in any longitudinally-Lorentz- 
transformed frame for p+p collisions in the PYTHIA 
model, consistent with the data from NA22 experiment. 
These results show that the string fragmentation mech- 
anism implemented in PYTHIA well describes the pro- 
duction mechanisms of charge particles and their charge 
balance in longitudinal phase space. 

Fig. 1(c) and (d) are the scaled balance function Bs{5rj) 
at two incident energies. They are deduced from di- 
rectly measured B[5ri\ri^) at six different pseudorapid- 
ity windows, (—0.8,0.8) (open circles), (1,3) (open tri- 
angles), (—3,1) (open squares), (—2.4,2.4) (open dia- 
monds), (0,3) (open crosses), and (-2,-1) (open stars). 
From the figures we can sec that all the Bs{St]) deduced 
from different windows are coincide with each other 
within errors, as expected from boost-invariance of the 
BF [4]. The solid down triangles in the same figures are 
the BF of the whole pseudorapidity range, B{6T]\rioo)- It 
is close to the scaled balance function Bs{6ri). These re- 
sults indicate that the scaled BF is in fact corresponding 
to the BF of the whole pseudorapidity range B{Sr]\oo) [3]. 

It is then interesting to see whether the boost- 
invariance of the BF is held in nucleus- nucleus collisions. 
STAR experiment only observe the boost-invariance of 



3 



9 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1 AuAu200GeV _ 

- AMPTvl.ll 

• (^3,-1) 

- C' ("2,0) — 

^ i A (0,2) 

- ' 1 ▲ (1,3) 

* 1 . 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1 AuAu 200GeV J 
- AMPTv2.11 

* ■ 

i 

-(b) 






J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 

- AuAu200GeV O (-0-8.0.8) I 

- AMPTvl.ll A Z 

1 1 V'^^'l 

(-2.4,2.4) - 
* (0,3) 

: (c) ® » > i i: 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4 

- AuAu 200GeV - 

- AMPTv2.11 - 

^ i ^ - 

T ^ P ^ 

" S B- 







0.5 1 1.5 2 0.5 1 1.5 2 



FIG. 2: Upper panel: the B{5rj\r]^x) in five pseudorapidity 
windows with equal size |77w| = 2 at the different positions 
for Au+Au collisions at = 200 GeV by (a) the AMPT 
default and (b) the AMPT with string melting. Lower panel; 
the scaled balance function, Bs{Sr]), deduced from the directly 
measured BF at various pseudorapidity windows with differ- 
ent sizes and positions for Au+Au collisions at y^s — 200 
GeV by (c) the AMPT default and (d) the AMPT with sting 
melting. 

BF in cental pseudorapidity range -1 < 77 < 1 [13], 
where the single particle distribution is ahnost flat, or 
boost-invariance. Now in model investigation, we can 
carefully examine the property in the whole pseudora- 
pidity range. 

The upper panel of Fig. 2 is the BF in five pseudorapid- 
ity windows with equal size 77^ = 2 at different positions 
(-3,-1), (-2,0), (-1,1), (0,2) and (1,3). Where the 
Fig. 2(a) and (b) are the results from the AMPT de- 
fault (vl.ll) and the AMPT with string melting (v2.11), 
respectively. Both figures show that the BF is boost- 
invariance in pseudorapidity range (-3, 3) in two versions 
of the AMPT. 

The lower panel of Fig. 2 is the scaled balance func- 
tions, which are obtained from directly measured BF at 
six different windows as indicated at legend of the figure, 
where the solid down triangles are the BF in pseudora- 
pidity range (-4, 4). It shows that the scaled BF does 
not depend on the size and position of the windows, and 
corresponds to the BF of the whole pseudorapidity in two 
versions of the AMPT, consistent with the results of p+p 
collisions in the PYTHIA model. 



charged particles under the constraint of global elec- 
tric charge balance. Global electric charge conservation 
not only applies to all final-state charged particles, but 
also constrains particles which are produced at the same 
proper time of evolution. It is argued that the transverse- 
momentum of final-state particles may be roughly used 
as a scale of the proper time of th eir p roduction in the ex- 
pansion of nuclear collisions [3, llSl . [l6l . [l3| ■ Examining 
the pt dependence of longitudinal property of the BF will 
provide direct access on whether particles in specified p^ 
range are consistent to be produced simultaneously with 
well balanced electric charge. 

So we turn to check whether the boost-invariant of BF 
holds for particles in different p^ ranges. Fig. 3 shows 
the BF for p+p colhsions at = 22 GeV and Vs = 200 
GeV from PYTHIA in three transverse momentum bins 
(0 < PT < 0.2), (0.2 <PT< 0.4), and (px > 0.2) GeV/c, 
respectively. These pt bins are selected to make the mul- 
tiplicity in each bin comparable. The result shows that 
the points at a given Srj in a restricted pt interval are 
approximately coincide with each other, i.e., the boost- 
invariance of the BF hold in small pt ranges. It indi- 
cates that particles produced at different pt ranges are 
also boost-invariant for hadron-hadron collisions in the 
PYTHIA model. 

The same study for Au-f Au 200 GeV collisions from 
the two versions of the AMPT are presented in the upper 
and lower panels of Fig. 4, respectively. Where four px 
bins are, (0.15,0.4), (0.4,0.7), (0.7,1) and (1,2) GeV/c. 
From the upper panel of the figure, we can see that the 
BF of different pseudorapidity windows in each pT bin 
are close to each other, in consistent with the data from 
STAR experiment However, in the AMPT with 

string melting, as shown in the lower panel of the figure, 
where the BF of different pseudorapidity windows are not 
as close to each other as those in the upper panel. 

This is because in the AMPT with string melting, each 
parton in the evolution of nuclear collision has its own 
freeze-out time, which last a very long period after the 
interaction of two nucleus . The particles in the same 
transverse-momentum range are not freezed-out simulta- 
neously with well balanced charge, and therefore the lon- 
gitudinal boost-invariance of the BF in small pT ranges is 
violated. In the AMPT default, the partons recombined 
with their parent strings immediately after they stop in- 
teracting, and converted to hadrons. So the charge bal- 
ance of the produced particles in the same ranges is 
preserved and boost-invariance of the BF keeps. 



V. SUMMARY 



IV. THE TRANSVERSE-MOMENTUM 
DEPENDENCE OF THE BOOST-INVARIANCE 
OF THE BF 

The longitudinal property of boost invariance of 
BF comes from the special longitudinal interaction of 



In the paper, we systematically study the longitudinal 
boost-invariance of charge balance function and its pt de- 
pendence for p+p and Au-|-Au collisions using PYTHIA 
the AMPT models. It shows that charge balance function 
is boost-invariance in both hadron-hadron and nuclear in- 
teractions, in contrary to the single particle density. As 
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FIG. 3: For each of three px ranges, the B{Sr]\Tjv,) in four pseudorapidity windows with equal size |77„| = 3 at the different 
positions for p+p collisions at ^/s — 22 GeV and ^/s — 200 GeV in upper and lower panels, respectively. 
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FIG. 4: For each of four px ranges, the B{5rj\ri„) in five pseudorapidity windows with equal size |t7w| = 2 at the different 
positions for Au+Au collisions at = 200 GeV from the AMPT default (in upper panel) and the AMPT with string melting 
(in lower panel) . 



expected, this boost-invariance of the BF make the BF 
properly scaled by window size is independent of window 
and corresponds to the BF of the whole (pseudo)rapidity 
range. Therefore, the BF is a good measure free from the 
restriction of finite longitudinal acceptance. 

It is further show that the boost invariancc of the BF in 
specified pt range is valid in PYTHIA for hadron-hadron 
collisions and the AMPT default for Au-fAu collisions. 
While the AMPT with string melting fails to reproduce 
this property due to the different schemes at hadroniza- 
tion. So the pt dependence of the longitudinal property 



of the BF may be served as a sensitive probe for charge 
balance in hadronization mechanism. 
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